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Abstract 

. 

C/3 ! The generic linear evolution of the density matrix of a system with a finite-dimensional 

qq ' state space is by stochastic maps which take a density matrix linearly into the set of density 

matrices. These dynamical stochastic maps form a linear convex set that may be viewed 
as supermatrices. The property of hermiticity of density matrices renders an associated 
' supermatrix hermitian and hence diagonalizable; but the positivity of the density matrix 

does not make this associated supermatrix positive. If it is positive, the map is called 
completely positive and they have a simple parametrization. This is extended to all 
Q> | positive (not completely positive) maps. A contraction of a norm-preserving map of the 

combined system can be contracted to obtain all dynamical maps. The reconstruction of 
the extended dynamics is given. 
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1 Introduction: Dynamical Maps 



The Quantum Density Matrix p is the statistical state and therefore quantum dynamics 
is the evolution of the density matrix. For a closed system this evolution is by a unitary 
time-dependent operator: 

p(t) = U(t)p(0)uHt). (1) 

The evolution is linear. But if we have an open system, its dynamics cannot be by a 
unitary evolution but a more general linear evolution [1]: 

p(t) = A(t) p(0) 

where A is a linear map of the density matrix into a density matrix. This superoperator 
A(t) can be written as a supermatrix. 

Pr,s(°) > A rS ,r's'(t)Pr's'(0) = (A(t)p) rs . (2) 

Then the supermatrix A(t) must satisfy the following three constraints: 

A S r,s'r' (t) ^rs,r's'{t) 

p(0) > — > A rs ^ s/ p r , s , > (3) 

-A-sss, r's' 8 r ' s f. 

These are consequences of the hermiticity, positivity and trace properties of the density 
matrices. These properties can be best expressed in terms of the dynamical matrix[2]. 

B rr i : s' s(t) = A- rs ^ r i s i(t) 

5 s's,rr-'(*) = B rr,s's(t) 

x*y r , B„>j s {t)x s y*, > 

B n r',s'n(t) = <W . (4) 
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Thus B is a hermitian matrix which gives nonnegative expectation value for super- 
vector that can be factorized: 

u^Bu>0 if u rs = x r y* . (5) 

It is sufficient if B is nonnegative but it is not necessary. If B > 0, then we will call the 
map "completely positive;" [2] but if only the positivity condition (4) is satisfied, we will 
call the map "positive but not completely positive" . Since B is hermitian, according to 
(4) it follows that it has an eigenvector decomposition 

B rr ',s's = £ Pa & = C (6) 

when M is the diagonal matrix with eigenvalue p a , and (^) are the normalized eigenvec- 
tors. For a completely positive map all the p a are nonnegative, but this will not be true 
for not completely positive maps. If all the p a are nonnegative, we can absorb them by 
defining the eigenvector 

r'O) - „ 1 /2 A<*) 

^rr' rex Vr' • 

So for, a completely positive map [2] 

p — ► Y C (Q) pC (a)t (7) 

a 

with the trace condition 

Y C^CW = 1 . (8) 

a 

2 Not Completely Positive Maps 

For a not completely positive map, some eigenvalues u are negative. We may then define 

p _^ C (Q)t P C (a) - Y DWp £>(/*) (9) 

a p 

with 

DfJ = (|,|) 1/2 ^ • (10) 



It will turn out that the number of negative eigenvalues cannot exceed the number of 
positive eigenvalue. The trace condition now becomes 

J2 - Yl Dm dW = 1 • ( n ) 

a (3 

C a u a ul C at - D®uf> u\ > . (12) 

a 

In dealing with maps we have the convexity property: 

B 1 cos 2 9 + B u sin 2 6 = B (13) 

is a positive map if B 1 and B n are positive maps. Every map that cannot be expressed 
in this manner is called extremal: 

B = B 1 cos 2 6 + B n sin 2 Q — > B 1 = B u . 

If we can determine all extremal maps, we can generate all maps from them. 

Since the density matrices are restricted by positivity and the trace condition, they 
form a compact set. The dynamical maps map from a compact set into a compact set, the 
maps also form a compact set.[l] Such a set can be generated from its extremal elements. 
Then 

]T tr (C^ C f ) - ]T tr (d^d^t) = N . (14) 

Define 

]T C (Q) C (a)t = J > , ]T £> {/3)t = K > . (15) 

a p 

We can do unitary transformations U, V on the initial and final density matrices we can 
diagonalize J. By virtue of (11) this will diagonalize K. Then 

J- K = 1 

If the eigenvalues of K are fc 2 and for J they are j 2 , then 

ka = 3a + 1 • 
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Define # a so that 



Define 



j a = sinhi? a , k a = cosh & a 



(16) 



C (a) = coshtf Q M(a) 
D {a) = smh$ a N(a) 



(17) 



so that 

m n 

E M W M U = 1 > E N(J3)N^(J3) = 1 • (18) 

a=l /3=1 

Parametrizising these matrices is already known [3]. In those cases where sinh^ a = 0, the 
matrix N(a) are not defined and we sum over a smaller set than the set of M{a). 

Choose a unitary transformation W\ so that is diagonal with eigenvalues 



cos 9\ 



(i). 



W\ M (1) Wi 



/ cos 0j 1} 



\ 



cos 0. 



(1) 



Than 



w\ c (a) C (Q)t Wi 



x 

/ sin® 0™ 



' • • cos 2 6$ ) 



sin 2 



V o 



\ 



'-.sin 2 ) 



Define 



M, 



(2) 
(a) 



/ sin 



V 



sm 



o(l) 



\ 



M(a), 2<a<N. 



•■smflW y 



Then define the (TV — 1) x (iV — 1) matrices which satisfies 



M^(m) M^m) = ljv-ixiV-i 



(19) 



(20) 



(21) 



We now repeat the procedure. Define 



W$M<®(2)W 2 



( cos 0f ] 



cos 8 



(2) 



\ 



• COS (fm ) 



(2) (2) (2) 

and a new orthogonal transformation a a such that #2 =0 an d we can define so 
that we can introduce the to — 2 matrices (m) which satisfy 



M< 2 >( 



/ cos o \ 



m = 



cos 4 3) 



(22) 



V ' • • COS® 3 / 

E M ( (3)m ) M £) = ^ 

This procedure can be carried out until all matrices are parametrized. The M x 
needs m parameter Q^p . . . , 9$ and so on until we get, a total of m 2 angles. The unitary 
matrices are also relevant parameters which determine the dynamical map. 

The same procedure can be carried out to parametrize K. We need n 2 angles. 
Together then, the matrices J, K together require angles d\ . . . , i? m , m 2 angles for J, n 2 
angles K; a total of (to 2 + n 2 ) parameters to determine the matrices C^ a \ to within 
to + n unitary matrices according to 

which leave (11) unchanged but change the 

3 Dynamical Maps as Contractions: 

A straightforward way [4] of generating positive maps is to consider a unitary evolution of 
a coupled system. R, S with S being the system and R a 'reservoir'. If we take a direct 
product density matrix 

p^pxr^VpxrV^ 
where V is a unitary matrix in the direct product space 



H m x H n . 
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Then 

Prs X Tab — > K'a.r'o' Pr's'Ta'b' V.*b,s>b ■ (23) 

The partial trace operation is a contraction 

V P X T -> y r „ jr ' a ' Pr's'Ta'b' V s b,s'b' (24) 

If r is made diagonal (if necessary by a unitary transformation in H n ) and the eigenvalues 
are r (1) , r (2) , . . . , r (n) , the map is 

Prs — > y rr'{n, v, r) p r ' s ' t{v) V* s , (n, v, t). (25) 

V 

r' s' 

n 

Clearly, if r has more than one nonzero eigenvalue, the map is not extremal. So as far as 
extremal maps are concerned 

Prs > V rr'{n) pr's' Kf s '( n ) ( 2 6) 

n 

which is of the standard form [3] 

p^J2 c (a) pC (Q)t 

a 

in which n a runs over at most 1 < n < m. We note that all these maps are completely 
positive. We can also do an inverse reconstruction: given an extremal completely positive 
map 

we can define a reservoir matrix r with all elements except t\\ zero and t\\ = 1. Then we 
can construct a unitary matrix V in mn dimensions with: 

Vra,r'l = C rr '(a) ■ (27) 

The conditions on (7(a) are transcribed into 

Y V ra,r', V s*a,s' > = $rs (28) 
a 
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which is necessary for V to be a unitary matrix. The ambiguity in constructing other 
elements does not affect the map. Thus a completely positive extremal map can always 
be obtained as a contraction of a unitary evolution[6]. The other elements can be 
used to carry out the construction of the matrix V and a generic diagonal matrix r. 

4 Not Completely Positive Maps as Contractions 

What about not completely positive maps? To obtain such a map by contraction we 
generalize the auxiliary space Hn to be a space with an indefinite metric and V to be 
a pseudounitary operator in the MN dimensional space. Positivity is guaranteed if the 
generalized density matrix is entirely within the convex of positive metric states of the 
MN dimensional space. Then the sum over n in (26) goes over both positive and negative 
metric terms; but the resultant density matrix is nonnegative. 

We can invert this derivation to realize the most general extremal not completely 
positive map as the contraction of a larger evolution in an indefinite metric space for the 
reservoir and a density matrix r of the reservoir system to have a single eigenvector (with 
positive metric) with eigenvalue unity, the others being trivial. 

Since such reservoirs are somewhat artificial, we may consider this reconstruction as 
a purely formal device. 
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Summary 



We have studied linear dynamical maps which take the set of density matrices 
into the set of density matrices. These maps form a convex set and are compact in 
the case of completely positive maps. The search for extremal maps gives us the 
restriction that we need at most N terms for an extremal map. Those maps can be 
obtained as contractions of a direct product system: the extremal maps correspond 
to a reservoir matrix which is a projection. Conversely we can reconstruct the 
unitary evolution of the expanded system from the map itself. 

The considerations are extended to positive but not completely positive dy- 
namical maps. The extremal maps still contain at most N terms. We can obtain 
these as contractions of an extended system with a pseudounitary evolution matrix. 
We could also reconstruct the extended pseudounitary evolution from the maps. 

In the systematic parametrization we find we need N(N — l)/2 parameters 
for the completely positive map apart from a set of unitary N x N matrices. For 
the corresponding parametrization of the not completely positive maps we have 
(m 2 + n 2 ) for m positive and n negative eigenvalues for the dynamical matrix and 
the unitary matrices (or less). 

These results generalize the results obtained two decades ago by Gorini and 
Sudarshan[6] for 2 x 2 matrices. 

Needless to say, however complicated the dynamical processes leading to the 
linear stochastic evolution that is represented by the dynamical map, we see that 
the same dynamics obtains with a reservoir having dimension iV 2 x iV 2 . In the case 
of an extremal map it suffices to have a reservoir with a state space of the same 
N x N dimensional density matrices. 

In this paper we have only dealt with dynamical maps, not the continuous 
semigroup of evolutions. This study was carried out by A. Kossakowski[9] and 
followed by others. [8,9] In these while the semigroup generators are parametrized 
no attempt is made to embed them in a larger system. Since the Zeno effect [10] 
operates for very small time intervals, care must be taken in generating a semigroup 
from the dynamics of an extended system. We hope to examine this question in the 
near future. 
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